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ABSTRACT 

This  paper  develops  the  mathematical  relations 
for  the  "Virtual  Kinematic  Model"  as  an  improvement 
over  the  vectorial  model  developed  earlier.  The  model 
solution  described  herein  provides  the  most  energy 
efficient  means  of  cooling  gears  -  that  is,  it  requires 
the  least  pressure  or  pumping  power  to  distribute  the 
coolant  on  the  tooth  surface.  Further,  this  nozzle 
orientation  allows  impingement  to  the  root  of  the 
tooth  if  needed  and  provides  the  most  cooling  control 
when  compared  to  into-mesh  or  out-of-mesh  cooling. 

NOMENCLATURE 

a  addendum  of  tooth 

dj  impingement  depth 

dx  uncorrected  impingement  depth  in  origin  plane 

ev  of  ©  =  sec  <t>  =  evolvent  function  of  9 

N  number  of  teeth  in  gear 

n  rotational  speed 

Pd  diametrical  pitch 

Pn  normal  diametria!  pitch 

Pt  transverse  diametrial  pitch 

R  standard  pitch  radius 

Rb  base  radius  of  gear 

R(jV  base  radius  of  virtual  involute 

Rj  radius  at  impingement  point 

Ror  operational  outside  radius  (for  modified 
addendum) 

R0  standard  outside  radius 


Rt  transverse  pitch  radius 

Rx  uncorrected  radius  at  impingement  point  (see 
Fig.  3) 

t  arbitrary  fl ight  time 

tf  jet  stream  flight  time  at  impingement 
tw  time  to  rotate  through  angle,  9W 
Vg  linear  pitch  line  velocity  of  gear 

V j  jet  velocity 

a  jet  incl ination  angle 

A a  addendum  modification 

d,  Pd 

Sj  dimensionless  impingement  depth  -■  - 

eiv  virtual  roll  angle  at  impingement  point 

eov  virtual  roll  angle  to  outside  diameter 

©a  vectorial  angle  of  involute  between  point  i 
and  0  (Fig.  2) 

9 j  involute  of  <j>i  =  involute  function  at  impinge 
ment  point 

9jv  virtual  vectoral  angle  to  impingement  point 
60v  virtual  vectora1  angle  to  outside  diameter 
0r  angle  across  tooth  space  at  base  circle 
0W  total  vectoral  angle  at  impingement  point 
v j  dimensionless  jet  velocity  Vj/Vg 

Pb  dimensionless  base  circle  radius 

p,  dimensionless  impingement  radius 


1 


o  dimensionless  outside  radius 

p  dimensionless  transverse  pitch  radius 

4>  pressure  angle  at  pitch  radius 

<ti  pressure  angle  at  impingement  point 

$iV  pressure  angle  of  virtual  involute  at  impinge¬ 

ment  point 

$n  normal  pressure  angle 

$0  pressure  angle  at  outside  radius 

<?ov  pressure  angle  of  virtual  involute  at  outside 

diameter 

transverse  pressure  angle 

41  helical  angle 

u  gear  angular  velocity 

g 


R i  after  time  (t  =  tf).  From  inspection  of  Fig.  1, 
it  can  be  seen  that  the  radius  of  impingement  is: 


Ri  -  R0  -  Vjtf  n> 

where  Vj  is  the  jet  velocity.  The  time  of  flight  of 
the  oil  jet  through  the  tooth  space  is: 


t 


f 


sec 


(2) 


The  gear  must  rotate  through  the  angle  0W  in 
the  same  time  it  takes  the  oil  jet  to  travel  tnrougn 
the  tooth  space  to  a  depth  of  dj  =  Vjtf  and  impinge 
at  Rj  .  From  Fig.  1  we  see  that  9W  =  ®o  +  ®r  +  -i- 
The  components  of  9W  are  developed  below. 

The  pressure  angle.  $i ,  at  the  depth  of  imoinge- 
ment,  di  is: 


u.  virtural  jet  angular  velocity 

jntrcdlcticn 


where 

Rq  is  the  gear  base  radius 


Cooling  of  lightweight  high  speed  gears  via  the 
lubricating  oil  supply  provided  to  the  gearbox  conti¬ 
nues  to  be  one  of  the  most  difficult  areas  in  the  high 
performance  mechanical  power  transmission  design  proc¬ 
ess.  Since  the  cooling  process  is  controlled  by  oil 
jet  orientation  and  by  the  jet  velocity  relative  to 
the  pitch  line  velocity  (PLV),  these  authors  have  con¬ 
sidered  the  development  of  a  fundamental  understanding 
of  the  oil  jet  primary  impingement  phenomena  as  essen¬ 
tial  to  understanding  how  the  gear  tooth  cooling 
process  works.  In  tne  early  1 970 ’ S  when  this  work  was 
started,  a  "vectorial  model"  was  used  to  describe  the 
radial  jet  flow  phenomena,  (Akin,  1975),  which  com¬ 
pared  reasonably  well  with  the  experimental  data, 

(Akin,  1975),  especially  at  pressures  above  approxi¬ 
mately  140  kPa  (20  psi).  Models  were  also  developed 
for  into  mesh  (Akin,  1983)  and  for  out  of  mesh  (Akin, 
1985;  to  determine  the  oil  jet  impingement  depths  for 
these  conditions.  The  into  mesh  and  out  of  mesh 
models  indicate  the  oil  jet  impingement  depths  were 
limited  and  therefore  maximum  cooling  could  not  be 
obtained  using  this  approach. 

The  "Virtual  Kinematic  Model"  was  first  mentioned 
by  Townsend  in  1980,  but  the  full  mathematical  devel¬ 
opment  of  this  model  has  not  been  published  until  new. 
This  rede1  has  been  generalized  to  include  spur  or 
helica!  gears  with  standard  AGMA  proportions  as  well 
as  long  or  short  addendums.  Further,  parametric  limi¬ 
tations  have  been  formulated  to  allow  the  gear  analyst 
to  know  he  is  specifying  his  problem  where  a  solution 
is  possible.  Also,  the  dimensionless  parameter  devel¬ 
opment  has  been  provided  to  allow  the  gear  engineer  to 
examine  a  wide  range  of  generalized  solutions  easily 
and  to  expose  the  limits  mentioned  above. 


Ri  =  R0  -  dj 

sc  that  6j  =  i nv  $j  =  tan  f  i  -  $>j  • 

Also,  the  angle  across  the  tooth  space  at  the  gear 
base  circle  is: 


9r  =  Jj-  -  2  inv  =  jj- -  2 C  tan  $  -  $)  (4; 

where  4>  =  cos'^R^/R)  is  the  pressure  angle  at  the 
standard  pitch  circle  R.  Further,  the  pressure  ang'e 


at  the  outside  radius 

Ro 

i  S : 

,  /R.'\ 

R  2 

*0 -  ccr  hr) 
\  0/ 

and 

tan  4>0  = 

-V  1 

[Rb 

so  that 

9  =  inv  d>  =  tan  *  -  4> 

0  0  0  0 

and  thus 

9  =  9  9  +  9.  as  stated  earlier,  so  that  tne  time 

w  0  r  1 

it  takes  the  gear  to  rotate  through  6  is 


t 


w 


30  9 

_ w 

rrn 


sec ; 


(5) 


where  n  =  rev/min. 


KINEMATIC  SOLUTION  FOR  MODIFIED  SPUR  GEARS 

The  jet  stream  flight  time  starts  (t  =  0)  at 
position  "1"  at  the  intersection  of  the  top  land  and 
trailing  profile  of  the  leading  tooth,  as  shown  in 
Ftg.  I.  Then  the  terminal  end  of  the  jet  passes 
through  the  tooth  space  where  it  impinges  at  "i"  on 
the  leading  profile  of  the  trailing  tooth  at  radius 


Equating  the  times  of  flight  t.  and  rotation  t  : 


d .  30  e 

1  _ w 

Vj  =  irn 


or 


Ro  -  Ri  !w 
vj  U9 


(7) 


If  dj  is  given,  then  the  required  V^  can  be  found 
expl  id  tly  as 


2 


Uqd.  ircLn 

0  =  30  6 

w  w 


note,  that  d.  is  a  function  of  $.(d.  =  f  (*.)) 

l  Mini 

If  Vj  is  given  and  tne  resulting  d^  or  R^ 
is  desired,  we  must  iterate  until  R.  ,  as  fn($j  ,  V. 
is  satisfied.  In  dimensionless  form,  where: 


=  Vj/Vg  and  5i  “  diPd/2 


vj  "  R 


oew  =  [ 


(M  +  2)6 


GEARS  WITH  MODIFIED  ADDENDUMS  AND  AN  OIL  JET  DIRECTED 
AT  AN  ARBITRARY  INCLINATION  ANGLE  a  FROM  THE  GEAR 
CENTER 

A  more  general  solution  can  be  provided  to  this 
problem  if  we  assume  that  tne  jet  is  pointed  away  from 
the  line  through  center  of  the  gear  by  an  arbitrary 
angle  "a"  as  shown  in  Fig.  2.  This  angle  is  defined 
to  be  positive  in  magnitude  when  pointed  in  the  direc¬ 
tion  of  wg  as  shown  in  Fig.  2.  Modified  addendum 
lengths  or  center  distances  are  accommodated  by  the 
use  of  the  parameter  aR0  =  ANg/2Pn  =  Aa,  so  that 
Ror  =  Rq  +  Aa.  Obviously  when  Aa  =0,  R0  =  Ror.  In 
this  paper  we  will  generalize  and  use  Aa ,  R0,  and 
Ror  in  all  equations . 

For  helical  gears  we  let:  Rj  =  N(2Pn  cos  <y) ,  so 


Ror  =  R0  +  Aa.  Obviously  when  Aa  =0,  R0  =  Ror. 


that,  R0  =  Rt  +  a 


=  R)-  +  a  +  Aa ,  where 


«  Pn  cos  +  and  the  subscript  t  refers  to  the 
transverse  plane  (normal  to  gear  axis.)  This  problem 
has  been  formulated  by  assuming  that  the  jet  terminal 
end  (or  head)  follows  a  trajectory  in  mesh  with  a 
''virtual  involute"  as  shown  in  Fig.  2.  The  virtual 
involute  is:  defined  by  the  radius  of  the  base  cir¬ 
cle  Rbv,  where 


■  Rcr  5in  «•  Vj 


“jRbV 


“j  =  Ror  5in  * 


As  can  be  seen  by  inspection  of  Fig.  2,  "u."  is 
the  theoretical  angular  velocity  of  the  virtual  jet 
involute  that  produces  the  proper  locus  for  the  term¬ 
inal  end  of  the  jet  stream  line  at  the  actual  jet 
velocity  Vj ,  subsequent  to  passing  the  point  1  at 
t  -  0.  The  virtual  involute  is  used  to  develop  the 
mathematical  relationships  needed  to  formulate  the 
time  of  flight  to  impingement  of  the  jet  stream  on  the 
trailing  tooth  profile  at  Rj  when  the  gear  angular 
velocity  is  ug.  Note  that  in  general  ug  *  u.  since 
the  virtual  base  circle  is  not  fixed  on  the  gear  blank 
(or  wheel ) . 

The  writer  has  selected  the  virtual  involute  gen¬ 
eralized  roll  angle  ev  as  the  parameter  to  interre¬ 
late  the  virtual  involute  for  the  jet  stream  with  the 
rotating  gear  geometry.  Thus  the  virtual  roll  angle 
eov  at  the  gear  outside  diameter  can  be  calculated 
from  the  fundamental  relationship  (see  Fig.  2): 

R,  c  =  R  COS  a,  SO  that  c 
bv  ov  or  ov 


From  here  we  can  develop  the  remainder  of  the  virtual 
involute  functions  including  60v  and  $ov  as  fol¬ 
lows  (see  Fig.  2): 

2a  ,,1/2  „  ,  2  , ,1/2 
e  =  (.ev  0-1)  or  ev  0  =  It  +1)  where  by 
ov  ov  3 

definition  (Vogel,  1945): 

ev  6  =  R  /R.  ,  so  that  c  =  [(R  /R.  )2  -  1]1/2 

or  bv  ’  ov  or  ov 

avoids  the  use  of  a  functions.  The  virtual 
vertoi  ial  angle  or  polar  angle  at  the  G.D.  is: 

6  =  [R  /R  )2  -  1]1/2  -  tan* 1 [ ( R  /R  )2  -  1]1/2 

ov  or  bv  or  bv 

(rad.)  and  the  virtual  pressure  angle  4>  is 
obtained  from: 

<j>  =  i nv* ' 6  =  sec"!(ev  6  )  =  cos*1  (R  /R  ). 

ov  ov  ov  bv  or 


From  the  above  we  cet  e  -6  +o  as  a  check. 

ov  ov  ov 

Also,  a  =  a  =  ir/2  -  4  (rad)  as  a  cneck) 
ov  Tov 

R  cos  a 

c  =  -a— — — -  =  COt  a  (0  <  a  <  f-  )  (15) 

ov  Rbv  2 

The  virtual  roll  angle  ejv  at  tne  impingement 
diameter  is  developed  as  follows: 


eO-i  =  eov  ”  civ  =  “j*f  =  R  sin  a 

o 

so  that  c,(tf)  =  Eov  -  Ujtf  =  (Ror  cos  a  -  v.tf)/ 
(Ror  sin  a)  =  cot  a  -  Vjtf/(RQr  sin  a). 

Note  from  Fig.  2  that  V^tf  =  d^/cos  a,  thus  when 
and/or  t^  are  unknown: 


R  cos  a  -  -  d. 

,  ,  ,  or  cos  a  _ ! _ 

Eiv  i  RQr  sin  a  a  RQr  sin  a  cos 


and  c  .  =  c  -  c • 
o-i  ov  iv 


cot  a  COt  a  +  R  sin  a  C0S  a  R  sin  2  a  <1j> 
Or  or 

Again  the  other  virtual  involute  functions  at  the 
impingement  diameter  are;  starting  with  the  vectoral 
angle: 


1  (rad. ) 


(R  -  d , ) 2 

1/2 

(R„„  -  d, )2 

or  i 

or  i 

R2 

Rbv 

-  id  n 

-  i 

bv 

R  COS  a 
Qr _ 


R„  COS  a 
or _ 

R„„  Sin  a 
or 


cot  a  (0  <  a  <  Y  > 


and  the  virtual  pressure  angle  at  Rf  : 


l(kl 


nr 


♦iv  =  C0S 


Further  noting  that 


< R  -  d.  > 
or  1 


R  -  d. 
or  1  , 


-  1  =  4>jv  +  6-v  as  a  check,  and 


However,  we  cannot  use  dj  to  calculate  0a 
except  as  an  approximation  because  the  angular  origin 
for  0a  is  along  the  radius  vector  Rx  (Fig,  3)  pass¬ 
ing  through  the  time  origin  at  "0"  while  dj  is  meas¬ 
ured  along  the  Rj  radius  vector  terminating  at  point 
"i"  on  the  leading  profile  of  the  trailing  tooth, 
which  is  the  terminating  vector  position  for  the 
angle  9a ;  thus  (see  Fig.  3): 

Rx  =  Rj  cos  6a  =  (Ror  -  dj)C0S  6a  (22) 


ai  v  "  f-  ♦iv(rad')  and  R.  -  Rbv(e;y  *  1)W2 


dx  -  Ror  -  Rx  =  Ror  ~  (Ror  -  dj)cos  9a  (See  Fig.  3) 
Noti"g  that: 


as  a  check.  Thus  the  time  of  flight  can  be  checked 
from : 


it  can  be  shown  that: 


tf  -  — 


RQr  sin  a  COS  at  pv  -  /  (  Ror  sin  a)J 


COS  a 


(sec)  (19) 


dx  =  Ror  C'2S  a  -  [(Ror  COS  a)2  - 

dj (2Ror  -  dj>]1/2  COS  a. 

so  that,  6a  =  tan-1  [dx  tan  a/(Ror  -  dx>j- 

Then  as  a  check  on  tne  calculated  value  of 


dj  -  ( Rq  *  oa)  -  Rj  =  Ror*  -  Rj 


d;  =  R„ 


R  -  d 
or _ x 


Treoreti cal iy ,  the  virtual  involute  is  a  dummy 
device  for  mathematically  describing  the  time  motion 
of  the  terminal  end  of  the  jet  stream  line.  As  such, 
tne  jet  stream  line  is  the  line  of  action  between  the 
virtual  involute  and  the  oil  jet.  Thus  the  coincident 
points  "i"  at  radius  rj  on  the  virtual  involute 
(v.i.)  and  Rj  on  the  gear  tooth  (g.t.)  define  the 
point  of  oi 1  impingement  at  time  tf  (v.i.)  =  tw 
(g.t.).  This  coincident  instant  in  time  is  the  result 
of  the  simultaneous  motions  of  the  (v.i.)  and  (g.t  ) 
from  time  t  =  0  when  the  (v.i.)  is  at  its  outer 
position  2  at  the  virtual  base  circle  (Fig.  2)  while 
tne  trailing  edge  of  the  leading  gear  tooth  is  at 
position  "O'1  where  tne  jet  line  crosses  the  O.D.  and 
the  (v.i.)  is  such  that  the  rotation  of  the  (v.i.)  at 
ug  without  the  relative  rotation  of  (u.  -  ug)  through 
time  tw  would  place  the  (v.i.)  origin  at  position  3 
at  the  virtual  case  circle.  Thus  the  (v.i.)  is  pre- 
ceived  to  rotate  faster  than  a.  by  the  amount 
(u.  -  ug>  such  that  it  rotates  at  ug  +  (u.  -  ug)  =  u. 
inJkeeping  with  the  jet  velocity  Vj  which  places  the 
(v.i.)  at  position  4  at  the  virtual  base  circle  and  at 
i  on  the  leading  profile  of  the  trailing  tooth  as 
Shewn  in  Fig.  2 

It  is  now  necessary  to  involve  the  gear  tooth 
geometry.  The  solution  will  involve  the  assumption 


r  i  =  Ri 


Ror  -  dj 


which  should  be  the  same  as  before. 

If  it  is  noted  that  R.  =  RQr 
(see  Fig.  2) . 


dj  then 


9.  =  inv 


v  <t> j  =  tan  <t> .  -  4>.  = 


-cos  1  rT  J 


Note  that  <j>.  is  a  function  of  Rj  and  dj  making  an 
explicit  solution  for  (Rj  or  dj)  impossible  when  Vj 
is  specified  leaving  dj  or  Rj  as  the  dependent 
variable.  Also. 


4^a  tan  * 

_  2.  _ _ Q 

0r  -  N  ~  R 


where:  9n  =  inv  $n  =  tan  <t>  -  for  spur  gears  and 

=  w/n  -  20a  tan  4^/R^  -  29^  for  helical  gears 


9.  =  inv  $  =  tan  $  -  <j>.  and  tan  $  =  tan  $  sec  w 

t  z  z  l  r  n 


and  Rt  =  2p“ 


■  Yf  C0S 


Further 


Thus  from  the  geometry  of  Fig.  2  it  can  be  seen  that 
9a  can  be  approximated  from: 


V.tf  sin  a 


d.  tan  a  (R„„  -  R.  )tan  a 
i  or  i _ 


a  "  Ror  -  Yf  co>  *  '  V  -  di 


9.  =  tan  <t>  - 

o  o 


Now  the  gear  angle  of  rotation,  9w,  between  points  1 
and  i  in  Fig.  2,  can  be  calculated  from 


9  =  9  *9  +  9  +  0 

w  o  r  i  a 


4 


so  that 


(38) 


3C(9_ 


e; 


t  =  —  = 
w  Ug 


ea) 


sec 


(3D 


where 


n  =  rev . /mi n 


It  is  now  possible  to  equate  the  time  of  flight 
tf  of  the  oil  jet  to  the  time  of  gear  rotation  tw,  so 
that 


306 


d; 


t  ( g  .  t . )  =  -  =  77; - r  =  tf( v  .  i . ) 

w  d  wn  (V.  cos  a)  f 


(32) 


Since  9.  is  fn(dj),  the  impingement  depth 

cannot  be  solved  explicitly  but  must  be  iterated 
numerically.  Thus,  we  solve  for  \T  so  that  an 

explicit  solution  is  possible  as  follows: 


V.  = 


ird  j  n 


d  ; 


i  306  COS  a  "  6  COS 
w  w 


(33) 


where 


9  =  6  +  6  +6,  +  9 

w  o  r  i  a 


(39) 


so  that  angular  geometry  expressed  in  fundamental 
terms  of  N,  AN,  4>,  41,  a  and  or  are; 

starting  with  the  vertorial  angle  8  (see  Fig.  1)  is: 

d 


tan  0^  =  d.  tan  5-  ,  so  that 
a  1  k  . 


6a  =  tan 

d 


-1 


45.  tan  a  COS  4/ 


N  +  2  COS  4/  +  2  AN  -  45(  COS  41 


Next  the  impingement  angle  6j  (see  Fig.  1)  is: 

/  2  \’/2 

-l[pb 


6.  =  tan  4.  -  =  I  -r  -  1  -  cos  — 


1  Ti 


V 


(41  ) 


It  would  be  desirable  to  graph  the  results  of 
these  formulas  as  in  Figs.  4  and  5.  Thus  the 
independent  variables  will  be  normalized  to  provide 
dimensionless  solutions  to  the  above  results.  The 
principle  parameters  used  are  Vj  and  dj  so  that 
we  will  express  them  in  terms  of  the  dimensionless 
parameters  ,  5)  and  N  as: 


v  -  .  A_ 

i  '  V  U  R 
J  g  g  or 


d. 


9  R  COS  a 
w  or 


and 


Vj  "  vj“gRor5i  "  WTOT 


A. 

2/P 


d;  P„ 


and 


(34) 


25. 


di  -  P- 


n 

(35) 


where 


thus 


and 


45  j  cos  41 


pi  =  1  '  N  +  2AN  +  2  cos  4- 
N  cos  4. 


Ph  = 


b  "  N  +  2AN  +  2  cos  4 


((N  t  2  COS  if  x  2AN)  -  45 ■  COS  4/) 
N  cos  4>. 


por  =  1 


(42) 


where  W.D.  is  the  whole  depth. 

;  time  of  fl ight  (t^ 

edge  of  the  leading  tooth  top  land  at  radius  R 
to  the  radius  of  impingement  Rj  on  the  leading 
profile  of  the  trailing  tooth  is: 

25; 


Thus  the  time  of  flight  (t^)  from  the  trailing 


or 


^f  V.OJ_? 


3  g  n  or 


cos  a 


(36) 


and  in  terms  of  N,  using 


(2Pn  cos  4) : 


t  e  - 


RQr  =  (N  +  2  cos  41  +  2 An ) / 


45.  cos  41 


f  =  [vjWg(N  +  2  COS  4/  +  2AN)C0S  aj 


(37) 


Now  this  is  valid  onl y  when  (N  is  given): 


Si  i 


N  +  2AN  +  2  COS  4;  -  N  cos  4> 


t 


4  COS  41 


or  when  5j  is  given 


N  >  2 


(25.  -  l>  cos  41  -  AN 


1  -  cos  4>. 


6. (max)  (43) 


=  N(min)  (44) 


the  tooth  space  angle  6r  at  the  base  diameter  is 
e  =  r-  nr  tan  -  29. 


r  "  N  N 


(45) 


where 


Then  the  time  of  rotation  t  from  the  trailing 

edge  of  the  leading  tooth  top  land  at  radius  RQr  to 

to  the  jet  line  when  it  intercepts  the  impingement 
po'nt  at  radius  RQr  on  the  leading  profile  of  the 

trai  1 ing  tooth  Is: 


9t  =  tan  4>t  -  4>t 


-  1 


1/2 


cos 


.1  ^b 
pt 


(46) 


where 


_ N _ 

pt  “  N  +  2  cos  y  +  24N 


and 


^  “  sec  *t 

and  finally  the  involute  angle  60  to  the  outside 
diameter  is: 


1/2 


0  =  i  nv  ±  =  tan  <t  -a 

o  cr  yor  yor 


-  COS  (pb> 


(47) 

again  it  is  possible  to  equate  the  times  of  flight  tf 
and  rotation  t  ,  now  in  terms  of  dimensionless  param- 

W 

eters,  so  that 


9 

w 


4  5.  cos  v 

Vj<N  +  2  COS  2.SN ) COS  a 


(43) 


Therefore,  if  the  dimensionless  depth  s(  <  Sj 

(max)  is  given  the  dimensionless  velocity  is  cal¬ 
culated  from: 


v 


j 


45.  cos  y 

6  (N  +  2  COS  t  +  2iN)C0S  a 

l  « 


(49) 


Also,  if  tne  dimensionless  velocity  v. 
the  dimensionless  depth  S*  is: 


is  given  then 


5; 


6wVj(N  +  2  COS  v(i  +  2<1N)  cos  a 
4  COS  '4/ 


(50) 


where  tne  limits  for  Sj  and/or  N  from  Eqs.  43  and 
44  must  be  observed. 

Further,  solving  for  Sj  must  be  done  by  an 
iteration  technique  since  9w($j>  is  a  function  of  Sj 
and  is  therefore  not  explicit. 


RESULTS  OF  COMPUTERIZED  PARAMETRIC  STUDIES 


As  noted  in  the  introduction  a  completely  new 
analytical  approach  has  been  taker,  to  development  of 
the  mathematical  model  for  what  we  identify  as  the 
< i nemat i c  mode  1 . 

The  effect  of  inclination  angle  (a)  in  the  Kine¬ 
matic  model  is  shown  in  Fig.  4  which  is  in  keeping 
with  the  physics  of  the  problem.  It  should  be  noted 
here  that  windage  has  been  neglected  in  keeping  with 
the  results  from  previous  studies 
(Akin,  1975). 

It  should  also  be  noted  that  when  fan  jets  are 
used  and  are  oriented  with  the  fan  perpendicular  or 
broadside  to  the  direction  of  flow  around  the  gear, 
the  resulting  impingement  depth  will  probably  compare 
better  with  a  curve  for  a  slighly  larger  inclination 
angle  than  indicated  by  the  pointing  direction  of  the 


jet  nozzle.  Also  the  effect  of  jet  velocity  on  depth 
Sj  is  shewn  in  Fig.  5.  The  curve  presented  here  is 
for  a  specific  gear  design:  pressure  angle  4>  =  20°, 
number  of  teeth  N  =  28,  inclination  angle  a  =  47° 
and  dimensionless  jet  velocity  Vj/Vg  =  v.  =  0.6463. 

Figure  6  (Akin,  1975)  shows  the  theoretical 
results  compared  with  the  experimental  data  gathered 
in  1974.  The  experimental  data  did  lie  substantially 
below  the  prediction  from  the  mathematical  models  a: 
low  oil  pressures.  It  is  assumed,  as  also  mentioned 
earlier,  that  this  is  the  result  of  (1)  the  windage 
effect  on  a  fan-jet  at  low  jet  velocity  and  (2)  pres¬ 
sure  loss  in  the  nozzle.  As  the  jet  velocity  is 
increased  the  experimental  data  approaches  the  curves. 

The  windage  study  conducted  in  Akin  (1975)  involved 
predicting  the  trajectory  of  the  oil  droplets  after 
they  passed  into  the  tooth  space  between  the  gear 
teeth.  Under  these  circumstances  tne  windage  affect 
is  considered  negligible  except  for  small  droplet 
sizes  around  0.001  cm  (0.0004  in.). 

CONCLUSION 

This  paper  develops  the  matnematica!  relations 
for  the  "Virtual  Kinematic  Model"  as  an  improvement 
over  the  vectorial  mede!  developed  earlier.  Tne  model 
solution  described  herein  provides  tie  most  energy 
efficient  means  of  cooling  gears-that  is,  it  requires 
the  least  pressure  or  pumping  power  to  distribute  tne 
coolant  on  the  tooth  surface.  Furthe-  this  nozzle 
orientation  allows  impingement  to  the  root  if  needed 
and  provides  the  most  control  when  compared  to  the 
into-mesh  or  out-of-mesh  cooling. 
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FIGURE  2.  -  KINEMATIC  MODEL  WITH  ARBITRARY  INCLINATION  ANGLE  "0". 
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FIGURE  3.  -  ENLARGED  DETAIL  OF  KINEMATIC  MODEL  TO  SHOW  GEOME 
TRY  CORRECTION  FOR  d,  AS  A  RESULT  OF  0,. 


FIGURE  4.  -  JET  ANGLE  a  WHICH  GIVES  MAXIMUM  DIMENSION- 
LESS  DEPTH  OF  IMPINGEMENT  6j  FOR  A  GIVEN  DIMENSION¬ 
LESS  JET  VELOCITY  Vj. 
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FIGURE  5.  -  DIMENSIONLESS  &j  FOR  A  GIVEN  DIMENSIONLESS 
JET  VELOCITY  Vj  (AT  OPTIMUM  VALUE  FOR  ALPHA)  a  =  47°, 
PRESSURE  ANGLE,  20°:  NUMBER  OF  TEETH,  28;  8  DIAMETRIAL 
PITCH. 
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FIGURE  6.  -  CALCULATED  AND  EXPERIMENTAL  IMPINGEMENT 
DEPTH  VERSUS  OIL  JET  PRESSURE  AT  4920  AND  2560  RPM. 
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